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ABSTRACT 

We  prove  an  existence  result  for  steady  flows  through  a  strip  which  can  be  regarded 
as  a  perturbation  of  rigid  motion.  The  fluid  is  viscoelastic  with  a  differential  constitutive 
law,  e.g.  an  upper  convected  Maxwell  model.  Particular  emphasis  is  focussed  on  the 
question  which  boundary  conditions  need  to  be  imposed  at  the  inflow  boundary  to  make 
the  problem  well-posed. 
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SIGNIFICANCE  AND  EXPLANATION 


\ 

Steady  flows  of  viscoelastic  fluids  can  not  be  uniquely  determined  by  imposing  bound¬ 
ary  conditions  only  for  the  velocities  as  in  the  Newtonian  case.  The  reason  for  this  is 
that  the  fluids  have  memory,  and  therefore  the  flow  inside  the  domain  is  affected  by  what 
happened  before  the  fluid  entered  the  domain.  This  leads  to  the  need  for  extra  boundary 
conditions  at  an  inflow  boundary.  The  nature  of  these  inflow  boundary  conditions  has 
not  been  analyzed  previously,  and  it.  is  certainly  dependent  on  the  constitutive  law.  In 
this  paper,  we  look  at  the  special  case  of  differentia)  constitutive  relations  with  a  single 
relaxation  mode.  We  consider  steady  transverse  flows  across  a  strip  which  are  small  per¬ 
turbations  of  a  flow  with  constant  velocity.  It  turns  out  that  in  this  case  two  extra  inflow 
boundary  conditions  are  required  in  two  dimensions,  and  four  in  three  dimensions.  This 
is  what  would  be  expected  from  an  analysis  of  characteristics,  but  it  contradicts  the  belief 
of  many  rheologists  that  it  is  possible  to  prescribe  the  extra  stress  at  an  inflow  boundary. 

The  problem  studied  here  is  of  potential  relevance  for  numerical  simulations  of  steady 
flows.  Many  of  the  flows  currently  simulated  are  on  infinite  domains.  Numerically,  these 
domains  are  truncated,  and  on  the  inflow  boundary  of  the  truncated  domain  people  usually 
prescribe  the  extra  stress.  According  to  the  analysis  in  this  paper,  this  is  an  overdetermined 
problem,  and  therefore  errors  must  be  expected  from  this  procedure  unless  the  artificial 
boundaries  are  chosen  far  enough  out.  Kc- Up-'Ct  c1  i'  -“ 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary  lies 
with  MRC,  and  not  with  the  author  of  this  report. 


INFLOW  BOUNDARY  CONDITIONS  FOR  STEADY  FLOWS 
OF  VISCOELASTIC  FLUIDS 
WITH  DIFFERENTIAL  CONSTITUTIVE  LAWS 


Michael  Renardy 


1.  Introduction 

Steady  flows  of  non-Newtonian  fluids  can  not  be  uniquely  determined  by  imposing 
boundary  conditions  only  for  the  velocities  as  in  the  Newtonian  case.  The  reason  for  this 
is  that  the  fluids  have  memory,  and  therefore  the  flow  inside  the  domain  under  consideration 
is  affected  by  the  fluid  motion  that  occurred  before  the  fluid  entered  the  domain.  This 
leads  to  the  need  for  extra  boundary  conditions  at  inflow  boundaries,  which  must  contain 
some  information  about  the  flow  history  outside  the  domain.  The  precise  nature  of  such 
inflow  boundary  conditions  is  not  understood;  it  is  certainly  dependent  on  the  constitutive 
law  of  the  fluid. 

In  this  paper,  we  consider  differential  constitutive  models,  which  relate  the  extra  stress 
tensor  to  the  velocity  gradient  by  an  evolution  equation  along  streamlines.  For  the  sake 
of  concreteness,  we  consider  the  special  case  of  an  upper  convected  Maxwell  fluid,  we 
emphasise,  however,  that  the  analysis  in  this  paper  does  not  depend  on  the  special  form  of 
the  nonlinearities  and  also  applies  to  other  differential  models  (for  examples,  see  (3], [5), [7]). 
In  a  steady  flow,  the  extra  stress  tensor  T  in  an  upper  convected  Maxwell  fluid  satisfies 
the  equation 


(u  •  V)T  -  (Vu)T  -  T(Vu)t  +  AT  =  qA(Vu  +  (Vu)T).  (1) 

Here  u  is  the  velocity  of  the  fluid,  and  r?  and  A  are  positive  constants.  The  Newtonian 
fluid  is  recovered  in  the  limiting  case  A  — *  oo.  This  equation  for  the  stress  is  to  be  solved 
in  conjunction  with  the  equation  of  motion 

p(u  •  V)u  =  div  T  -  Vp  +  /,  (2) 

div  u  =  0. 

The  fact  that  (1)  describes  an  evolution  of  the  stress  along  streamlines  may  suggest 
that  appropriate  boundary  conditions  for  (1)  and  (2)  are  given  by  prescribing  T  at  an 
inflow  boundary,  in  addition  to  the  Dirichlet  conditions  for  the  velocity  which  are  required 
in  the  Newtonian  case.  In  fact  this  seems  to  be  what  rheologists  generally  believe  (see 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041  and 
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e.g.  j2],  p.  31).  The  study  of  characteristics  (4),|6),j9),  however,  suggests  otherwise.  We 
may  regard  (1)  and  (2)  as  a  first  order  system  for  p  and  the  components  of  u  and  T. 
In  this  system,  it  turns  out  that  the  streamlines  of  the  flow  are  double  characteristics  in 
two  dimensions  and  quadruple  characteristics  in  three  dimensions,  while  the  number  of 
independent  components  of  the  symmetric  tensor  T  is  three  and  six,  respectively.  From 
this  it  appears  that  prescribing  the  stress  at  an  inflow  boundary  may  overdetermine  the 
problem. 

In  an  earlier  paper  [&),  the  author  has  proved  an  existence  result  for  slow  steady 
flows  of  differential  non-Newtonian  fluids.  These  flows  were  in  a  bounded  domain  with 
prescribed  homogeneous  Dirichlet  conditions  on  the  boundary.  In  this  case,  there  are  no 
inflow  boundaries,  and  no  extra  boundary  conditions  are  needed.  The  existence  proof 
was  based  on  an  iteration  scheme  that  alternates  between  solving  an  elliptic  system  and 
a  hyperbolic  system  which  has  the  streamlines  as  characteristics  (a  similar  approach  was 
recently  used  by  Beirao  da  Veiga  [lj  to  prove  the  existence  of  steady  flows  of  compressible 
Newtonian  fluids). 

In  the  present  paper,  we  consider  transverse  flows  through  a  strip  which  are  small 
perturbations  of  a  flow  with  constant  velocity  and  zero  stress.  If  we  want  to  use  the 
iteration  scheme  of  [8j,  we  have  to  impose  boundary  conditions  at  inflow  in  order  to  be 
able  to  solve  the  hyperbolic  part.  If  this  is  done,  it  turns  out  that  in  fact  the  iteration 
converges,  but  in  general  not  to  a  solution  of  the  equations  (1)  and  (2).  An  extra  condition 
has  to  be  satisfied  for  this  to  be  the  case.  We  shall  then  modify  the  iteration  scheme  to 
accommodate  this  condition.  If  this  is  done,  the  number  of  inflow  conditions  required  is 
as  the  analysis  of  characteristics  suggests. 


2.  An  iterative  procedure  for  steady  flows  perturbing  flow  with  constant 
velocity 


We  consider  transverse  flow  through  the  strip  0  <  x  <  1  with  periodic  boundary 
conditions  imposed  in  the  y-  and  2-directions.  The  periods  are  denoted  by  Ly  and  Lz. 
For  /  -  0,  a  solution  of  (1)  and  (2)  is  given  by  a  flow  with  constant  velocity  U  >  0: 
u  -  (V.  0, 0),  p  ~  0,  T  =  0.  We  are  looking  for  flows  which  are  a  small  perturbations 
of  such  a  flow.  I.e.,  we  prescribe  a  small  body  force  /  and  boundary  conditions  for  the 
velocity  which  are  small  perturbations  of  the  constant  velocity:  u(0,y,2)  =  (l/, 0,0)  + 
Vi(y, z),  u(l,y,s)  -  (f/,0,0)  +  v2(y,2).  These  boundary  conditions  must  be  consistent 
with  incompressibility, 


*z  •  v\(y,  z)  dy  dz  = 


ez  •  v2{y,z)  dy  dz. 


(3) 


In  addition,  we  shall  need  conditions  on  the  stresses  at  the  inflow  boundary  x  -  0.  The 
nature  of  these  conditions  will  now  be  discussed. 

Following  |8],  we  apply  the  divergence  operator  to  the  constitutive  equation  (1),  and 
obtain 

(u  •  V)div  T  -  (Vu)div  T  *  Adiv  T  -  T  :  d2u  t)XAu.  (4) 
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Here  we  have  set  T  :  d2  =  £ .  k  j-yd,:a-yt  ■  If  we  substitute  div  T  from  the  equation  of  motion, 
we  find 

V|(u  •  V)p  +  Ap]  -  IVu  +  (Vu)TjVp  -  (u  ■  V)/  +  (Vu)/  -  Xf 

=  T  :  (?2u  +  qAAu  -  p(u  •  V)(u  ■  V)u  -r  p(Vu)(u  •  V)u  —  Xp(u  ■  V)u.  (5) 

In  [8],  the  construction  of  a  solution  was  based  on  an  iterative  method,  which  alternates 
between  solving  a  “Stokes-like”  problem  and  a  hyperbolic  equation  whose  characteristics 
are  streamlines.  Similar  iterations  are  in  fact  used  in  numerical  calculations  (2j.  Since  we 
are  looking  for  small  perturbations  of  the  flow  with  constant  velocity,  we  can  take  this 
flow  as  a  starting  point  for  the  iteration.  The  iteration  scheme  of  |8]  is  the  following  (we 
introduce  the  new  variable  q  =  (u  ■  V)p  +  Xp) 

u°  =  (f/,  0,0),  p°  =  9°  =  0,  T°=  0,  (6a) 

T"  :  dV+1  -  qXAun+1  -  p(u"  •  V)(un  •  V)tin+1  -  Ap(un  ■  V)un+1  -  V9"+1 
=  -j(Vu")  -i-  (Vun)T]Vpn  -  (un  •  V)/  +  (Vun)/  -  Xf  -  p(Vun)(u"  •  V)un, 
div  untl  =  0.  un+)  =  (17,0,0)  +  U)  on  x  =  0,  un+1  =  (17,0,0)  t  iij  on  i  =  1, 

f'fL*[L'qn+'dzdydx  =  0,  (66) 

Jo  Jo  Jo 

(u”+1  •  V)p”+1  +  Apn+I  =  9n+1,  (6c) 

(t£n+1  •  V)Tn+1  -  (Vun+1)Tn+1  -  Tn+1(Vu’>+,)T  +  ATn+1 

=  rjAKVu'1-”1)  +  (Vun+1)Tj.  (6d) 

In  setting  up  this  iteration,  it  is  essential  that  the  singularly  perturbed  operator  (u  ■  V)  +  A 
is  inverted  rather  than  evaluated.  As  long  as  U*  <  Xrj/p,  and  un  and  T"  are  small 
perturbations  of  (U,  0,0)  and  0,  equation  (6b)  is  an  elliptic  system  for  un+I  and  qn+1 .  If 
l/1  >  A  p/p,  a  change  of  type  occurs  |4],[6],[10],  and  this  case  will  not  be  considered  here. 
Equations  (6c)  and  (6d)  are  hyperbolic  equations  with  the  streamlines  as  characteristics.  If 
there  are  no  in-  and  outflow  boundaries  as  in  [8],  characteristics  do  not  cross  the  boundary, 
and  (6c)  and  (6d)  can  be  solved  without  any  boundary  conditions  imposed.  In  the  present 
situation,  however,  we  have  to  prescribe  p  and  T  at  the  inflow  boundary  x  =  0  in  order  to 
solve  (6c)  and  (6d). 

It  can  in  fact  be  shown  along  lines  similar  to  |8)  that  with  such  inflow  boundary  condi¬ 
tions  the  iteration  will  converge,  provided  that  the  data  are  sufficiently  small.  However,  it 
will  in  general  not  converge  to  a  solution  of  the  original  problem.  The  “solution”  obtained 
from  the  iteration  will  satisfy  (1)  and  (5),  but  the  original  problem  is  (1)  and  (2).  Let  us 
recall  that  (5)  was  obtained  by  substituting  p(u  ■  V)u  +  Vp  -  /  =  div  T  into  (4),  and  (4) 
was  obtained  by  differentiating  (1).  Hence  the  fact  that  (5)  is  satisfied  means  the  following: 
If  we  set 

p(u  ■  V)u  +  Vp  -  /  -  div  T  =  h,  (7) 

then 

(u  •  V)h  -  (Vu)fc  -  Xh  =  0.  (8) 
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If  streamlines  do  not  cross  the  boundary  as  in  |8j,  then  (8)  does  in  fact  lead  to  h  =  0, 
but  in  the  present  case  we  need  to  know  that  h  =  0  at  the  inflow  boundary  in  order  to 
make  that  conclusion.  This  condition  has  to  be  viewed  as  a  constraint  on  the  possible 
inflow  data  for  p  and  T,  and  the  iteration  (6)  has  to  be  modified  in  order  to  accommodate 
this  constraint.  If  we  simply  count  numbers,  we  see  that  in  three  dimensions  T  has 
six  components,  p  has  one  and  h  has  three,  thus  suggesting  that  four  inflow  boundary 
conditions  can  be  prescribed.  In  two  dimensions,  T  has  three  components,  p  has  one 
and  h  has  two,  suggesting  that  one  needs  two  inflow  boundary  conditions.  This  is  also 
suggested  by  the  analysis  of  characteristics  [4),[6j,j9j. 

We  want  to  restrict  the  choice  of  inflow  boundary  conditions  in  (6c), (6d)  by  the 
constraint 

div  T""1  -  Vpn+1  =p(un+1  •  V)un+1  -  /  at  i  =  0.  (9) 

On  the  left  hand  side  of  (9),  the  ^-derivatives  of  pn+1  and  Tn+I  can  be  expressed  using 
(6c)  and  (6d),  e.g. 

dpn+l  A  „»+i 

— - —  =  -  —  pB+1  +  - 1-  nonlinear  terms.  (10) 

ox  U  U 

When  doing  this,  (9)  assumes  the  form 

A(pn+>  _  rr,+1)  +  T£yl  +  r,t;  = ... 

-  jj 1  +  (r2v‘  -  P”-')y  +  Tg*  = ...  (11) 

.  tjtn-fl  I  trpn+l  n+ 1  \  

fj1  13  1  23, jf  "r  U  33  ~P  )*  ~ 

The  right  hand  sides  indicated  by  dots  contain  terms  involving  /,  un+I,  qn'l'i  and  nonlinear 
terms  which  also  involve  Tn+1  and  pn+I. 

We  have  to  use  equation  (11)  in  order  to  express  some  stress  components  at  x  =  0 
in  terms  of  others.  To  obtain  a  convergent  iteration,  we  want  to  do  this  in  such  a  way 
that  no  loss  of  regularity  occurs  when  solving  for  the  undetermined  stress  components. 
Unfortunately,  this  is  not  possible  in  such  a  way  that  certain  components  of  stress  are 
prescribed  at  inflow  and  others  are  left  to  be  determined.  Rather,  we  have  to  prescribe 
stress  components  partially.  For  example,  we  can  do  the  following.  Let  each  component 
of  the  stress  be  expanded  in  a  Fourier  series  in  y  and  z ,  e.g. 

Tu(0,y,z)  =  (12) 

k,l 

We  can  then  prescribe  the  Fourier  components  of  Ttl,  T22,  Tl3  and  T33  for  |fc|  >  |f|,  k  f  0 
and  solve  for  those  of  p,  T\i  and  T33.  For  |/|  >  lk\,  we  prescribe  the  Fourier  components 
of  r,  it  7*121  r22  and  J33  and  solve  for  those  of  p,  T13  and  733.  For  k  ~  l  =  0,  we  prescribe 
the  Fourier  components  of  Tu,  7*22,  T23  and  T33  and  solve  for  those  of  p,  7*|2  and  Ti3.  In 
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the  two-dimensional  case,  this  amounts  to  prescribing  7n  and  and  solving  for  p  and 

r,  2. 

The  solution  procedure  is  now  as  follows.  We  start  the  iteration  with  initial  data 
(6a).  Then  at  each  step  of  the  iteration  we  compute  a  new  u  and  q  using  (6b).  Then  we 
compute  a  new  p  and  T  from  (6c)  and  (6d)  with  inflow  boundary  conditions  which  are  in 
part  prescribed  and  in  part  computed  from  (11)  in  the  manner  outlined  above.  We  show  in 
the  next  section  that  such  an  iteration  converges  under  appropriate  smallness  conditions 
for  the  body  force  /,  the  velocity  boundary  data  i/j  and  vj  and  the  prescribed  part  of  the 
stress  boundary  data. 


3.  Proof  of  convergence 

As  usual,  we  denote  by  H‘  the  space  of  all  functions  on  the  strip  0  <  x  <  1,  which  are 
periodic  with  periods  Ly  and  L,  in  the  y-  and  ^-directions  and  have  s  derivatives  which 
are  square  integrable  over  one  period.  Sobolev  spaces  of  periodic  functions  living  on  one  of 
the  boundaries  x  —  0  or  x  =  1  are  denoted  by  The  corresponding  norms  are  denoted 

by  ||.||,  and  ||.||(#). 

In  the  following,  let  s  be  any  integer  >  1.  We  assume  that  the  body  force  and  the 
velocity  boundary  data  satisfy  the  bound 

ll/IU+l  <  Tf>  ll»lll<.+3/2>  ^  11^2 1| (#+3/2)  <  It  (13) 

where  -y  is  a  positive  number  which  will  later  be  chosen  small.  The  stress  at  the  inflow 
boundary  consists  of  a  prescribed  part  Tp  and  an  unknown  part  Tu,  which  must  be 
determined  from  (11)  at  each  stage  of  the  iteration.  We  assume  that 

l|Tp||<.+,>  <  r  (14) 

The  convergence  proof  consists  of  two  parts:  First  we  show  that  all  iterates  remain 
bounded  in  a  certain  norm,  and  then  we  use  this  fact  to  prove  convergence  in  a  weaker 
norm.  This  procedure  is  typical  in  dealing  with  hyperbolic  problems,  and  the  iteration  (6) 
involves  the  hyperbolic  part  (6c), (6d).  In  order  to  carry  out  the  first  part,  let  us  assume 
that 

||u"  -  (f/, 0,0)||, +2  <  £,  llpnSU+.  <  €,  IIT1.+  1  <  £,  (15) 

where  c  is  small.  As  long  as  pU2  <  r>X  and  c  is  sufficiently  small,  (6b)  is  an  elliptic  system 
for  un+1  and  qn+1,  and  a  standard  argument  shows  that  there  is  a  constant  C\  such  that 

ll«n+I  -  (l/, 0,0) II. +2  +  ||vn+1||»+i  <  Ci (7  +  e2)  =:  6.  (16) 

In  the  next  step,  X”+ 1  and  pn+1  on  the  inflow  boundary  are  determined  from  (11).  Using 
the  trace  theorem,  we  see  that  qn+l  and  Vun+1  are  in  By  using  this  in  (11),  it 

is  easy  to  show  that  for  small  6  there  is  a  constant  C j  such  that 

I|T:+,||(.+1>  +  l|pn+lH(.+  t)  <  C2(7  +  6)  =:  o.  (17) 
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We  now  turn  to  the  hyperbolic  equations  (6c)  and  (6d).  The  solution  to  these  equations 
can  be  obtained  by  integrating  along  streamlines  (since  t/n+1  is  Lipschitz,  streamlines  exist 
and  are  unique).  To  obtain  estimates  for  the  solution,  let  us  multiply  equation  (6c)  by  p 
and  integrate  over  the  domain.  This  yields 


/  /  ’  f  '  pnJri(un~i  V)p”  +  'i-r\(pn+')2dzdydx  =  /'  [L=  p"+V+I  dzdydx. 

Jo  Jo  Jo  Jo  Jo  Jo 

(18) 

The  left  hand  side  is  equal  to 

^  J  /,  ** '  u"+,(l’V’j:)(Pn^ J',(l,!/,«))2  dz  dy 

f  f  ez  ■  un*l(0,y,z){pn~l(0,y,z))2  dz  dy  +  \  [  f  [  (pn  +  1)2  dz  dy  dx. 

*  Jo  Jo  Jo  Jo  Jo 

(19) 

The  first  term  in  this  is  positive,  and  hence  we  obtain 

A  f  f  f  (pn+1)2  dz  dy  dx 
Jo  Jo  Jo 

-JJ  J(i  dz  dy  dx+  l-  J  J  ez  •  un+1(0,y,z)(pn+1(0,t/,z))2  dz  dy. 

. J"  °  °  (20) 

That  is,  we  have  an  estimate  of  the  L2-norm  for  p  in  terms  of  the  i2-norm  of  q  and  the 
L2~norm  of  the  inflow  boundary  data.  By  differentiating  (6c)  with  respect  to  y  or  z ,  we  can 
obtain  estimates  for  derivatives  of  p  in  these  directions,  and  estimates  for  i-derivatives  can 
be  obtained  from  the  equation  itself.  Similarly,  we  can  deal  with  (6d).  If  6  is  sufficiently 
small,  this  argument  leads  to  estimates  of  the  form 


||pn+l!l.+  ,  +  ||T"+'|l.+  .<C3(a  +  ',  +  «).  (21) 


If  we  now  choose  t  small  enough  and  ~i  sufficiently  small  relative  to  t,  we  will  have  6  <  e 
and  C3(o  +  T(  +  6)  <  e.  This  implies  that  (15)  holds  with  n  replaced  by  n  +  1.  By  induction, 
we  see  that  (15)  holds  for  all  values  of  n. 

The  fact  that  all  the  iterates  are  bounded  and  in  fact  small  can  now  be  used  to  show 
convergence  of  the  iteration  in  a  weaker  norm.  More  specifically,  we  can  show  that  u" 
converges  in  H,+  l  and  that  pn,  qn,  Tn  converge  in  He.  The  argument  is  fairly  routine, 
and  we  shall  just  show  a  typical  step.  Let  us  take  (6c)  at  step  r  +  1  and  step  n  and  take 
the  difference.  We  find 


(un+1  ■  V)(pn^  p")  +  A(p"+1  -  p")  +((«n+l  -  un)-V)pn  =  qn  +  l  -  qn.  (22) 
From  this  we  find 


ilPn+,-Pn!l.  <C(!i,nJ|  +  V 


u  !ls+iitP  i's+i  +  \p 


i<«> 


).  (23) 


6 


T 


From  above,  we  already  have  a  bound  on  the  term  ||pn||»+i-  Similar  arguments  are  applied 
to  (6a),  (6d)  and  (11).  By  putting  the  resulting  estimates  together,  we  obtain  an  inequality 
of  the  form 


||u"+1  -  u"L+i  +  ||qn  + 1  -  qnt  +  |!Pn+1  ~  Pn|l*  +  l|Tn+I  -  Tn\i, 

<  R'(t,7)[||«n  -  +  ',9*  -  q”~ 1 IU  ^  IIP"  -  p"_1ii<  +  l!T"  -  T"-'|U].  (24) 

The  constant  K(e,  7 )  tends  to  zero  as  7  and  1  tend  to  zero.  Hence  a  contraction  is  obtained 
if  we  choose  e  and  7  sufficiently  small. 
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